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Predictability and semiclassical approximation at the onset of 

black hole formation 



We combine analytical and numerical techniques to study the collapse of 
conformally coupled massless scalar fields in semiclassical 2D dilaton grav- 
ity, with emphasis on solutions just below criticality when a black hole almost 
forms. We study classical information and quantum correlations. We show ex- 
plicitly how recovery of information encoded in the classical initial data from 
the outgoing classical radiation becomes more difficult as criticality is ap- 
proached. The outgoing quantum radiation consists of a positive-energy fiux, 
which is essentially the standard Hawking radiation, followed by a negative- 
energy flux which ensures energy conservation and guarantees unitary evolu- 
tion through strong correlations with the positive-energy Hawking radiation. 
As one reaches the critical solution there is a breakdown of unitarity. We show 
that this breakdown of predictability is intimately related to a breakdown of 
the semiclassical approximation. 
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I. INTRODUCTION 



Quantum radiation from black holes is necessary in order to maintain the consistency 
of the second law of thermodynamics with the existence of black holes 0. On the other hand, 
evaporation of the black hole reveals one of the most fundamental problems in theoretical 
physics; the question of unitary evolution of the Universe. Does evolution from an initial 
pure state take place non-unitarily to a final mixed state [Q], or unitarily to a final pure 
state 0] ? One of the major obstacles to a better understanding of the Hawking effect is the 
complexity of four-dimensional (4D) semi-classical gravity f^. Simplified models which may 
give insight into the possible answers are two-dimensional (2D) dilaton quantum gravity 
theories 0. The dilaton field, viewed as part of the geometrical structure restores dynamics 
in 2D analogous to that of spherically symmetric 4D Einstein gravity. We consider the 
formation and evaporation of a 2D black hole by the collapse of massless matter scalar 
fields. The evaporation of the black hole via production of quanta of the matter fields can 
be fully traced in the 2D semiclassical theory, including the back-reaction of the evaporation 
on the geometry [[7|-[TT[]. 



If the energy and energy density of the infalling matter are sufficiently large, then the 
incoming matter forms a black hole. Otherwise, the original incoming matter escapes to 
infinity and no black hole is formed. In this latter case the evolution is unitary and no 
information is lost. These unitary solutions of the semiclassical theory are called subcritical 
solutions [|T^-|TB[. The study of subcritical solutions just below the critical threshold in 
which a black hole is formed may help us to understand the process of semiclassical black 
hole formation and its infiuence on information. Moreover, the interesting results obtained in 
classical gravity concerning critical behavior at the onset of classical black hole formation | T6[ 



make it important to examine this critical behavior in the context of semiclassical physics. 
In the critical regime of the 4D Schwarzschild black hole, the black hole mass approaches 
zero and the curvature near the horizon becomes large, so semiclassical effects must be 
considered. Since the 4D semiclassical theory is quite complicated, insight may be gained 
by considering a 2D theory that shares many of its dynamical features, namely, 2D dilaton 
gravity . The 2D dilaton models can be derived from 4D almost extremal dilatonic black 
holes [0 using the Kaluza-Klein reduction . 



Although the general properties of the 2D subcritical solutions, namely that they are 
stable and unitary, have been known for sometime [p!^ -[l^, investigation of the explicit 



evolution reveals some new features of physical importance |]T5[ . In Ref . [T^ we studied the 



evolution of the subcritical solutions with infalling matter in the form of shock-waves. In 
this work we extend our study to include the behavior of subcritical solutions for general 
smooth initial data, with emphasis on the near-critical solutions. 

In Sec. II we present our model of 2D semiclassical dilaton gravity as an initial value 
problem. We also derive the general equations to be integrated numerically for arbitrary 
initial data. In Sec. Ill we give examples involving smooth infalling matter. The previous 
shock-wave results |]T5| appear as a limiting case of these examples. 

In Sec. IV we address the question of information. The information that may be lost 
in the process of black hole evaporation is related to the correlations between the outside 
world and the interior of the black hole. Two types of information are involved: (i) "Classical 
information", carried by the classical matter that forms the black hole, and (ii) "Quantum 
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information", encoded in quantum correlations between outgoing and incoming pairs of 
particles created by the collapse geometry. A quantity that plays an important role in 
understanding the structure of the subcritical solutions is the outgoing radiation reaching 
future asymptotic null infinity, 3"^. In previous work an explicit form of that radiation 



was found for the first time in semiclassical gravity. This outgoing radiation is intimately 
related to the question of information. In principle, for the subcritical solutions one should 
be able to recover the complete information given by the initial data from that outgoing 
radiation. We show how in the subcritical solutions, classical and quantum information is 
encoded in the outgoing radiation reaching asymptotic future null infinity. 

In Sec. V we consider the subcritical solutions just below criticality. We show that as in 
the classical case, also in the semiclassical case solution space can be divided continuously 
into two regions, i.e., there exist continuously varying parameters, pi, in solution space, such 
that for Pi < p* the evolved scalar field will not form a black hole (the subcritical solutions), 
while for pi > p* a black hole will be formed (the supercritical solutions). We show that 
as the critical solution is approached (p, p*) the outgoing energy flux diverges and the 
fluctuations in the outgoing energy density become very large, implying a break down of the 
semiclassical approximation at criticality. 

In Sec. VI we show that near criticality the density of information encoded in the out- 
going radiation reaching becomes very large and diverges at criticality. This divergence 
results in an apparent breakdown of predictability that coincide with the breakdown of the 
semiclassical approximation. We present our conclusions in Sec. VII. 

II. THE MODEL 
A. 1-loop effective action 



Recently we have proposed a modified theory of 2D semiclassical dilaton gravity pO 
The effective action of the modified theory is 

SqQ = ScGHS + NSpL + Scott, (1) 
where Scghs is the Callan-Giddings-Harvey-Strominger (CGHS) classical action |T^, 

N 



1 „ 1 " 

Scghs = ^J d\^g ^e"^^ (p^^ + 4(V0)^ + 4A^) - \ ^(V/, 



(2) 



SpL is the Polyakov-Liouville action pO[ that incorporates the 1-loop corrections correspond- 
ing to the trace anomaly of the stress energy momentum tensor of each of the N quantum 
matter fields, 

n 

and 



SpL = -— (fx^-g{x) / (fx'^-g{x')R'^^\x)G{x,x')R'^^\x'), (3) 



/ d'xV^ {{V<l>r - (4) 
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is a local counter-term that we add in order to get an exactly solvable theory. In the 
above is the dilaton field, i?*-^-* is the 2D Ricci scalar, A is a positive constant, V is the 
covariant derivative, and G{x, x') is an appropriate Green's function for V^. The real value 
functions, fi{x), are the classical values of the massless scalar fields. One can regard each 
of the fi{x) as the expectation value of the quantum field operator, fi{x), in an appropriate 
quasi classical coherent state, |a) The effective action (|l|) describes the full quantum 



theory in the large N limit, in which case the fluctuations of (f) and g^^, can be neglected 
||l^ , p!0[| . Recently Mikovic showed that one can derive the effective action of Eq. (|I]) from 
ScGHS by flxing the diffeomorphism gauge, solving the constraints, and then quantizing 
the reduced system. After choosing an appropriate initial quantum state, one recovers the 
action (P l-loop effective action |22| . 



In null coordinates, 2;^, and conformal gauge, (7++ = g = 0, (7+- = — (1/2) exp(2p), 

the action (|l]) takes the form 



S'gg = — I dz^dz 



(9_F)9+(X - -Y) + (9+F)9_(X - -Y) + A2exp(-2r) 



2 . V ^ / V 2 



1 ^ 



i=l 



+ (5) 



where X = exp(— 20), Y = 0— p, and k, = Nh/12. In the large N limit we take h to approach 
zero while keeping k flnite. The kinetic action density of the system described in Eq. (|5|) is 
a bilinear symmetric form (9+G)M(0)(9_G), where G is a vector comprised of the {N + 2) 
flelds X, Y and the N matter flelds, /j, and M{(f)) is an {N + 2) x {N + 2) symmetric matrix. 
One can verify that the determinant of M is proportional to = exp(— 40), and unlike in 
other models of modifled dilaton gravity 0, here this determinant is non- vanishing for all 
real values of 0. The vanishing of the determinant at X{x^, x~) = signals a singularity. 

The equations of motion derived from varying the action (|^) with respect to X, Y and 
fi are 

d+d-X = -A^ exp(-2y), (6) 
d+d.Y = 0, (7) 
d+d.f, = 0. (8) 

The constraints (from varying the action ([l|) with respect to g±±) are 

- dlX - 2d±Xd±Y - T^^ + k [{d±Yf + dlY + t±{z^)\ = 0, (9) 

where after varying with respect to gf^^ {fi = ±, u = ±) we set g++ = g = to get 

(I). Here T^^. = {l/2)J2i{d±fiY is the classical contribution to the energy-momentum 
tensor of the N matter flelds, and tj-{z^) are integration functions determined by the 
speciflc quantum state of the matter scalar flelds. A conformal coordinate transforma- 
tion of the form z^ — > and z~ —>■ y~ , preserves the form of the metric, g±± = and 
= — (1/2) exp(2p), where the new conformal mode function is related to the old one 
by p{y^, y^) = p{z^, z^) + In \J dz^ /dy+ + In ^dz~ ldy~ . The dilaton fleld is a scalar, and 
therefore the fleld Y transforms like —p. The general solution of Eq. (0) for Y is Y{z^, z~) = 
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Y^{z^) +Y_(z^), and we can choose the coordinates = exp[—2Y±{z^)]dz^ , for which 
Y{x^, x^) = 0. In the following we use these "Kruskal" coordinates, denoted by (x"*", for 
which (f){x^,x~) = p{x~^,x~). Henceforth, the indices ± on tensors will refer to components 
in the Kruskal coordinates. 

In the Kruskal gauge the general solutions to (|^) and (|[), subject to the constraints 

are 



X(x~^,x ) = —X'^x'^x — 



dxt T^\{xi)-K,t+ 



dx^ Tf'_ (x^ ) — Kt_ {x^ ) 



(10) 



and 



/,(a;+,x-) = /+(x-) + /r(x-) 



(11) 



B. The initial value problem 



Consider first the linear dilaton (LD) solution which corresponds to fi{x^,x~) = and 
t±(a;^) = in (p!OD, namely, the solution Xld{x^, x~) = —X^x^x~ . It is defined for < < 
cxo and —oo < x~ < 0. In the manifestly flat coordinates, cr^ = r ±cr = ±A~^ ln(±Ax^), the 
LD solution corresponds to the fiat metric, ds"^ = —da~^da~, and the dilaton field has the 
linear form, (p = —Xa. As shown in Fig. |I|, the null curve x"*" = defines left asymptotic past 
null infinity, the null curve x~ = —oo defines right asymptotic past infinity, The null 
curve x~ = is left asymptotic future infinity, and x"*" = +oo is right asymptotic future 
infinity, ^3^. In general, the initial data on and determine completely the solution 
in the region < x"*" < oo, — oo < x~ < 0. Specifying these initial data is equivalent to 
giving /j"(x^), and t_(x^) on and fj^{x~^) and t+(x"'") on Giving these functions, we 
integrate (|T0]) to find the solution everywhere, however such a solution may not be physically 
acceptable in the whole space-time, since singularities may appear. One can say that our 
1-loop effective theory is exactly solvable however there is one major difficulty with such an 
approach. Consider first the LD solution: while on we have a oo and exp(20) = 0, 
on ^2 ~^ — oo) we have exp(20) oo. From the action (|^) we see that exp(20) plays 
the role of the "coupling constant", and so the coupling diverges on S^^. One can split the 
space-time into a region of weak coupling and a region of strong coupling, see Fig. Those 
two regions are divided by a curve, the "boundary curve" (specified below in more detail). 
In the strong coupling region we cannot trust the 1-loop effective theory, especially on 53^. 
Therefore giving the initial data on is in general unphysical. One way to avoid this 
problem is to consider the solutions only in the weak coupling region and impose on a time- 
like curve boundary conditions that preserve unitarity and conserve energy. These criteria 
are satisfied by imposing reflecting boundary conditions on the boundary curve. Thus, the 
initial value problem that we define is the following: on we give the initial data, f^{x~^) 
and t+(x"'"), and on the boundary curve defined by ) = p(x ), we impose the 

reflecting boundary condition. 
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Before specifying the boundary condition we elaborate on the energy-momentum tensor. 
The c- number T^^, which one gets after varying the effective action in (|l|) with respect to 

the metric g^'^, can be regarded as the expectation value of the quantum operator, T^^, in 
the quasi-classical coherent state, \a), Namely, 



where in the expansion of {a\Tj[j^\a) the first term, T^l, is of order and is the classical 
part of the stress tensor, while the second term, (T^j^u) , is of order h (or A^^ for fileds) and 
is the one-loop contribution. The one-loop contribution to the energy-momentum tensor is 
given by 

(T,±,±) = 4dl^p - {d,±pf - t,±(z±)] ; (T,+,-) = Kd,+ d,-p , (12) 

in some general null coordinates {z^, z~). Equation (|12D follows by integrating V^(T^'^) = 
and using the trace anomaly (T^^) = for the massless matter fields fi |^ 



or equivalently by varying the Polyakov-Liouville action (^ with respect to the metric 

One can argue that operationally the split between the classical and quantum radiation 
is not well-defined when we consider a single (scalar) field. However, for fields one could 
classically excite only some of the fields. Then the total radiation in the remaining fields 
is just the quantum part. To operationally distinguish among the A^ fields, one can, for 
example, add another quantum number to these scalar fields. 
The reflecting boundary condition is |25jT3|j5|,p!5[] 



tL{x-) = {p'{x-)fTl+{xUx-)) + K{p\x~)f/^dl{p'{x-))-^/^ (13) 

where ' = d/dx~ . The last expression on the right hand side (r.h.s.) of Eq. (|13]) is due to 
quantum particle creation from the boundary, which is effectively a moving mirror . 
One can split the reflecting boundary condition (|T^ into its classical and 1-loop parts, 

T^'_(x-) = (p'(x-))^T;V(4(x-)) (14) 

Ux-) = {p'{x-))H4xUx-)) + T/i^; 15 

If we take the flelds fi to satisfy Neumann or Dirichlet boundary conditions, then the classical 
reflecting boundary condition is satisfled. 

Next we consider the boundary curve. We would like it to separate the regions of weak 
coupling , exp(20) < g"^, and strong coupling, exp(20) > g"^, where g"^ is some parameter that 
specifles the value of the coupling below which the 1-loop effective theory is trustworthy. 
Since X{x^,x~) = exp(— 20), we can deflne the boundary curve to be the curve on which 
X{x~^, x~) = Xb = g^"^ = const. The solution X{x^ , x~) should be determined by the initial 
data on and by the boundary condition (|1^). The boundary curve, X(x+,x~) = X^, 
depends on the solution X(x'^,x~). Thus, we get a highly non-linear problem, unlike the 
straightforward problem, with initial data given on ^3", having solutions (|T0|) and (p!T|). We 
next reduce this non-linear problem to solving a single ordinary differential equation. 
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C. The boundary equation as a second order ODE 



The function T^_^_{x~^) can be viewed as the part of the initial data on that describes 
the classical profile of the infalling matter. We take T^'+(a;^) to be a general function of 

with compact support. On the other hand, t^{x~^) describes the quantum state on $5^. 
Since in this work we would like to study the Hawking effect, we take the quantum state 
to be such that on 53];^ we have no quantum radiation, (i.e., {T^y){v) = 0, where v is the 
asymptotically flat null coordinate on Q]^, v = ln(Ax"'")). Since in the asymptotically 



flat null coordinates the conformal mode, p, approaches zero on Q"^, we get from ([T^) that 
in these coordinates ty{v) = 0. To find the corresponding t^{x^) in Kruskal coordinates, 
we use the tensor transformation of {Tz±z±) in Eq. (^2|) (under a conformal coordinate 
transformation) and get 



dx+ ; V ' ' 2 V (2x+) 
where -D^f-z] is the Schwarz operator defined as 




D'yiz] = {diz)/{dyZ) - ^-[d'yZ/dyZ] (l7) 

and we use t„(f) = 0. Eq. (|l^) then becomes 

{p'{x-)y/'dUp'{x-))-'/' 

+ 1 + n/{2Xs) • ^''^ 

The second term on the r.h.s. of Eq. (p!8| ) is the result of particle creation from the dynamical 
boundary. Recall that in the 4D Hawking effect m , the "boundary curve" is r = (the fixed 
point of spherical symmetry). The curve r = does not act as a moving mirror and there is 
no particle creation at r = 0. The creation of particles (i.e., the Hawking radiation) is due 
to the curvature of space-time near the horizon. Can we also eliminate the moving mirror 
effect in our 2D theory? We see from (|18D that if Xb << k, then the moving mirror term 
is negligible. This is consistent with the fact that X = exp(— 20) is indeed the 4D radial 
coordinate |T8|p6| , and Xb corresponds to r — in 4D. However, one should be careful 



when taking the limit Xb ^ 0. It corresponds to the limit g'^ —>■ oo, which defines the strong 
coupling region. The 1-loop effective theory (|1]) is trustworthy as long as g^h << k << 1. 
In this case (in the region of interest, i.e., X{x'^,x~) > Xb) the quantum corrections for 
the dilaton-gravity part are negligible compare to the 1-loop corrections for the N scalar 
fields, which in turns are small compared to the classical contribution. Therefore, we need 
to satisfy both g'^h « k, « 1 and Xb « n- With = g^, we combine the conditions 
and get 

4 « Xb « Nh « 1. (19) 

One can always take the large N limit in such a way that (0) is satisfied and Xb is arbitrarily 
close to zero. In the following we take this large N limit. Eq. (|18|) then reduces to 
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1 / p'{x 
4 \ p{x~ 



(20) 



One can write the general solution (10) in the form 



X{x+, X-) = -Xx+iXx- + A-^P+(x+)) - ^ ln(Ax+) + ^^^^ + F{x-), 

4 A 



(21) 



where M{x~^) and P^{x~^) are the mass and momentum of the infalling classical matter, 



M(x+) = X xtT^+ixt)dxt 

J 



PAx^ 



(22) 



and 



Fix- 



doc n 



cia;]^ T'^__{xi ) — Kt_{xi 



(23) 



Here F(x~) is a function of x~ to be determined by the boundary conditions. Before 
we obtain the equation for the boundary curve, let us define the following dimensionless 
quantities: 



z = \x~ 



and 



q{z) = \x%{x (z)) = \p{x (z)) 



(24) 



where q{z) is a dimensionless function that specifies the location of the boundary curve. 
Using (^), (pOD, (pl|) and acting with d'^/dz'^ on the boundary equation X{x%,x~) = Xb, 
we get a second order ODE for q{z) 



z + Pg{q) + 



K 



4g(z) 



d'^q ^dq ^ 
dz^ dz 



4:q^{z) 



dz 



0, 



where 



Tf^iq) ^ \'T++ixUx-)) and P,(g) = f Tf^{q)dq 



(25) 



(26) 



are given functions of g, which are determined by the initial data on SJ^. Eq. (|25|) is a 
second order, non-linear, ordinary differential equation for the boundary curve q{z). We 
solve it numerically for different profiles of the infalling matter using an embedded fifth- 
order Runge-Kutta ODE integration routine After obtaining q{z) we use Eqs. (|T^ and 
(13) to find Tf'_(x~) and t-{x~) and integrate Eq. (p!0|) to obtain the solution X{x^,x~). 



III. SOLUTIONS 



A. Vacuum solutions 



In this subsection we consider solutions of ( P^D with different profiles of infalling null 
matter. First consider the vacuum solutions, for which T^Hq) = Pq{q) = 0. The general 
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solution of 

{dq/dz){z - 



with T^j^iq) 



Pq{q) = and the initial conditions q{z — oo) — > and 



— oo 



0, is 



where a is some positive constant. One finds from 
of vacuum solutions 



(27) 

2]) and (^7|) the one-parameter family 



— ln(— Ax^x ) + — ln(a) 



a. 



(28) 



These are static solutions |T^. For a > k/A the solution (|28| ) has a time-like singularity in 
the strong coupling region, e"^^ —>■ oo, while for a < k/A the solution has null singularities in 
the weak coupling region e"^"^ — > 0. The solution with a = k/4 is everywhere regular, has the 
geometry of a semi-infinite throat, and can be regarded as the ground state of the theory 
|lO| . One can show that compared to this ground state, the ADM mass of the static 
solution in Eq. (ESl) is 



M = A 



K , /4a 
— In — 

4 \k 



K 

"-4 



(29) 



B. Smooth infalling matter 

In the following we consider more general infalling matter with compact support, Xi < 
x~^ < X2- For x+ < xf, i.e., region I in Fig. |^, we have T^^ = Pg = 0, and the solution is one 
of the vacuum solutions (^). If the total mass of the infalling matter, M = M{x~^ oo) 
in Eq. (p2D, is above a critical value Mcr, then a black hole is formed. The critical mass 



Mcr is of the order of Xa and is determined numerically from the profile of the infalling 



matter stress-tensor and the ADM mass of the initial spacetime. We take M/ X^ k for the 



validity of the semiclassical approximation |29|. Since we consider sub critical solutions in 



which a black hole does not quite form, a is at least of the order of M/ X and hence k. 
We take a = 1 and k = 10~^. Our results would be similar for any values of a and k as long 
as a ^ K. We also fix the scale of the x^ coordinate such that Xxi = 1. Define zi to be 
the value of z = Ax" in Eq. (pSj) corresponding to q{zi) = Xxf (see Fig. j^). Using Ax^ = 1, 
a = 1 and Xb ~ 0^ we find from Eq. (p8|) that zi = —1. For z < zi, i.e, in regions I, II 
and III in Fig. 0, the solution of Eq. (^) for q{z) is given by Eq. (|27D, and the solution for 
X(x+, x~) is 

X(Ax- < zi) = -Xx+[Xx- + A-^P+(x+)] - y ln(-A2x+x~) + ^^""^^ 



A ' ' X 



For z > Zi = —1 (see Fig. ^ we have to integrate Eq. ( P^D numerically to find the boundary 
curve q{z). So the range of z to be integrated numerically is [—1,0). The initial values for 
the numerical integration of the ODE (pSj) are (using (|27|) with a = 1): q{z = zi) = 1 and 
dq/dz{z = zi) = 1. 
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Before we give the results of our numerical integration, it is convenient to define coordi- 
nates in which the dynamical metric of our spacetime corresponding to the solution (^Tj) is 
manifestly asymptotically flat on Q"^. In the Kruskal coordinates, the metric is 

ds^ = — exp[2p{x'^ ,x^)]dx~^dx~ = —X~^ dx~^dx~ , (30) 

where X = exp[—2p{x~^,x~)] in these coordinates and is given by Eq. (pi|). The required 



asymptotically fiat coordinates are (m, v), defined in terms of the Kruskal coordinates by the 
following conformal coordinate transformation: 

V = X-^ \n{Xx+) and u = -X'^ ln(-Ax" - P+/A) , (31) 

where P+ = P+(x^ oo) is the total momentum of the infalling matter. Rewriting the 
metric (^) in terms of the {u, v) coordinates shows that exp[2p(f , u)] — ^> 1 as — > oo on 
and u — >• —oo on $5^. 

Consider now the stress tensor describing a smooth profile of infalling matter: 

, V <0 

T^liv) = { (2M/e) sm^i-Kv/e) , < i; < e (32) 
, f > e , 



where v = ln(Aa;"'") . In Fig. ^ we show the profile (|3^). Here M is the total mass of 
the infalling matter, M = T^l{v)dv, and e is its width. 

To integrate ( pSj ) we need to write T^^ in terms of q. Using ( P^ ) and ( P^ ) in (P^), we get 

r , ln(g) < 

T^M = ^sin2(vrln(g)/(Ae)) , < ln(g) < Ae (33) 
( , ln(g) > Ae . 

We numerically solve (^) with ( |33D for different values of M / A and Ae, and find the corre- 
sponding boundary curves q{z) and solutions X{x~^,x~). 

It is hard to extract physical information from the coordinate-dependent definition of the 
boundary curve. For example, as long as the boundary curve is everywhere timelike, one can 
define null coordinates = x"*" and y~ = x^{x~) in which the boundary curve is "static", 
being given by the equation y = 0, where y is the spatial coordinate y = {y^ — y^)/2. On the 
other hand a quantity giving physical insight into the nature of the solutions is the outgoing 
radiation reaching future asymptotic null infinity, The outgoing radiation consists of 
a classical part ([T^) that is refiected from the boundary, and a quantum part ([T^). We 



calculate both in the manifestly asymptotically fiat null coordinates on 53^. In terms of the 
u coordinate defined in Eq. (|3l|) we have on Q"^, 

T:l{u) = iXx- + P^/XrTUx~{u)) . (34) 

On the other hand, Eq. (|12D gives 

{Tuu){u) = -Ktu{u) , (35) 
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where we used the fact that p and its derivatives vanish on in the manifestly asymptoti 



cally flat coordinates. Furthermore, using the coordinate transformation (0) and Eq. (16) 
with — s> x~ and v ^ u, one can express tu{u) in terms of t_(x~(-u)) to get 

{T^u){u) = (/€AV4)[1 - \-\\x- + P+/\)H_{x-{u))] . (36) 

Putting the numerically integrated function q{z) (defined in Eq. (0)) and the initial data, 
T^_^{x~^), in (0) and (pOD, we find the profiles of the outgoing matter fluxes T^^ and (Tuu)- 
The results are shown in Fig. ^ (a) and (b) in which the infalling matter is described by the 
profile given in Eq. (|32|) with Ae = 5. In Fig. ^ (a), the mass is M/X = 0.1 (which is 10% of 
a), while in Fig. ^ (b) it is M/X = 1. The solid curves depict the total stress-tensor of the 
outgoing radiation, = T^l + {Tuu), in units of A^, while the gray curves correspond to 
the quantum radiation stress-tensor, (Tuu)- Since k, = 10~^ is much smaller than M/X, the 
quantum radiation is scaled up in the figures by a factor of 2/k = 2 x 10^ and hence does 
not show up in the solid curve. 

As we will see later in section V, when the total mass and the local density of the infalling 
matter are large enough, then a black hole is formed; otherwise the infalling matter escapes 
to infinity without forming a black hole. For a fixed width e in Eq. (|32|), the mass of the 
infalling matter M fully determines its profile. Let M^r be the mass above which a black 
hole is formed. As discussed in section V, for Ae = 5, we find that M^r/X ~ 1.5. In Fig. 
§ (a), we take the mass of the infalling matter to be about 6% of the critical mass M^r- 
We see that the total outgoing radiation is very similar to the incoming radiation. Fig. |^a. 
The profile of the matter field is almost unaffected by the evolution. The curvature and 
coupling are everywhere small and the infalling matter is just refiected from the boundary 



with almost no distortion. In the 4D case, this corresponds to "weak" initial data for 
which the null (scalar) field goes through the origin of the radial coordinate, r = 0, without 
developing large densities, and escapes to infinity with very little distortion. 

The quantum part of the outgoing radiation stress-tensor (T„„) in Fig. ^ is almost sym- 
metric. (This is true only for a symmetric incoming classical profile far from the critical 
point.) The quantum radiation includes a region of negative-energy density that insures en- 
ergy conservation. In section IV we will calculate the quantum correlation function and show 
that the negative-energy radiation is strongly correlated with the positive-energy radiation. 
This strong correlation is necessary for the quantum state on to be pure. 

In Fig. ^ (b) the mass of the infalling matter is about 65% of the critical mass Mcr- 
We see that the outgoing classical radiation is more distorted compared to that in Fig. ^ 
(a). In this case the infalling matter strongly distorts the spacetime before getting refiected 
from the boundary and escaping to infinity without forming a black hole. Likewise the 
quantum radiation in Fig.|^(b) is no longer symmetric. The reason for this is the following: 
As M increases, the total amount of positive-energy quantum radiation increases. Because 
of energy conservation also the total amount of negative-energy radiation increases with 
M. While the width of the positive-energy radiation increases with M, the width of the 
negative-energy radiation. Am, satisfies the quantum inequality [|T^,0 

\Eneg\^U < K , (37) 

where Emg is the total amount of negative-energy quantum radiation and Am is its width 
(as measured by an asymptotic observer). Since \Eneg\ increases with M, the width Am must 
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decrease and we get a non-symmetric profile of (Tuu)- As we will see in section V, as M 
approaches Mcr the width Am approaches zero, and we get a brief burst of negative energy. 

As a consequence of general covariance of the effective action, Eq. (0), we have V^T^^, = 
0. We find numerically that the energy defined in terms of Tj[^ is conserved, i.e., /q- T^^dv = 

Jcj+ Tl^du. This is what one would expect if the reflecting boundary at X = exp(— 20) = 0"*" 
corresponds to a static boundary like r = in 4D, as discussed in the paragraph after Eq. 
(p!8|). The h expansion implies that also V^T^J, = and V^{T^^) = (the latter was used 
in arriving at (T^,^) from the trace anomaly). Nevertheless, we find numerically that the 
energies associated with the classical and quantum parts alone are not separately conserved. 
For example, Jc^- T^idv ^ T^^du. Evidently, because the total energy-momentum tensor, 

R R 

T/.,, determines 0, and hence the boundary, it behaves like a static boundary only for 



/if 



but not for T^l or {T^u) individually. Far below criticality, the classical and quantum parts of 
the energy are nearly conserved separately. However, as one approaches the critical solution 
the classical and quantum parts of the energy each are strongly non-conserved. 
Next we consider the profile described by 



M 
ivo+e/2) 





sin^(7rf /e) 

M 



M 
{vo+e/2) 



vo+t/2 



sm' 



[ti{v - \vQ)/e\ 




t; < 
< i; < e/2 

I < t; < f + I 
fo + f<t;<fo + e 
f > f + e . 



(38) 



In Fig. |^(b) we show the profile (^Sj) with v^/e = 100, which is nearly a square waveform. 
For Vq ^ e, as shown in Fig. |5(b), the energy-momentum tensor (^) describes infalling 
matter of nearly homogeneous density. The total width of the profile is [vq + e). In terms 
of the g-coordinate defined in Eq. (0), the profile (^) is 



M 





sin2[7rln(g)/(Ae) 



M 



M 

A2(t;o+e/2)g2 



A2(i,o+e/2)g2 

sin^{7r[ln(g) - 




Xvo]/{Xe)} 



ln(g) < 
< ln(g) < Ae/2 
Ae/2 < ln(g) < X{vo + f ) 
X{vo + f ) < ln(g) < X{vo + e) 
ln(g) > X{vo + e) . 



(39) 



We take (pOf ) with Xvq = 10 and Ae = 0.1, and integrate (pSf ) for different values of M/X. 
Using the numerically integrated function, q{z), we calculate the outgoing classical and 
quantum radiation on shown in Fig. ^ (c) and (d). In (c) we take M/X = 0.1 and 
in (d) M/X = 2. Because the mass distribution of (^) is more spread out than that 
of (|3^) , Mcr is larger, having the value M^r/X ~ 3.0. In these figures, the solid curves 
describe the total radiation while the gray curves describe the quantum radiation scaled 
up by 2 X 10^. In Fig. ^ the mass is about 3% of the critical mass, which is the case of 
weak initial data, and is similar to Fig. The total outgoing radiation is very similar to 
the incoming classical radiation. Fig. ^(b), and the quantum radiation is spread across the 
outgoing classical radiation symmetrically. In Fig. ^ the mass is about 65% of the critical 
mass, as was the case in Fig. ^d. The outgoing radiation is more distorted and the quantum 
radiation is non-symmetric. 
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Finally, recall that in Ref. |T5 



the case in which the classical matter had the form of an 

consider profiles 
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incoming shock-wave was studied. To recover the results of Ref. 
with M/A = 0.1 but with different values of the width Ae. In Fig. ^, the curves a and a 
are, respectively, the total and quantum outgoing radiation in the case Ae = 5. The curves 
h and [3 are the total and quantum outgoing radiation for Ae = 1, and c and 7 correspond 
to Ae = 0.1. One can see that as Ae becomes smaller, the (incoming and outgoing) classical 
radiation becomes very localized, while the transition of the quantum radiation from positive 
to negative-energy values becomes more sudden. In the limit Ae — > 0, the classical energy- 
momentum tensor can be described by a delta function T^'^ = {M / \xq)5{x'^ — x^) and 
one recovers the results of Ref. ||T5[. One should remember however that as Ae 0, the 
derivative p' in Eq. ([T5|) diverges, and the large limit should be taken in such a way that 



the second term on the right-hand-side of Eq. (115]) can be neglected |]T5 



IV. INFORMATION 

A. Classical Structure 

In section III we considered infalling matter with very little structure. In order to encode 
non-trivial information in the classical infalling matter one should consider more complicated 
profiles. For example one can encode information in bits. In this case the infalling null matter 
can be a sequence of pulses each corresponding to one bit of information. Let us send a 
message in binary numbers. The number "1" is described by a pulse of relative height 1 
and "0" by a pulse of relative height 0.5. As an example consider the following incoming 
energy-momentum tensor 

16 

Tt{v) = Y.kQ,{v), (40) 



i=l 



where 



1 , elsewhere , ' 



k G {0.5, 1, 0.5, 0.5, 1, 0.5, 0.5, 0.5, 0.5, 1, 1, 0.5, 1, 0.5, 0.5, 1}. (42) 



and 



In Fig. ^c we show the profile (^Ol). It corresponds to the binary number 0100100001101001, 
which is the word "Hi" in ASCI. In the g-coordinate the profile (|40D is 

16 



r.1b] = E^^0^(?)' (43) 



where 

( 2OM 



Q.iqs = J UM^ sin {107r[ln(g) - i]} , i< ln(g) < z + 0.1 
* 1 , elsewhere . 
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Using (^31), we integrate (pSj) for different values of M/X and find q{z). Then we calculate 
the classical and quantum radiation on Q"^, shown in Fig. In (a) we take M/A = 0.1, and 
in (b) we take M/X = 2. The solid curves describe the total outgoing radiation, while the 
gray curves describe the quantum radiation, scaled by 2/k = 2 x 10^. 

For the profile (|iO|), we find numerically that the critical mass is Mcr/X ~ 2.5. In Fig. 
Ha the mass is about 4% of the critical mass, and indeed we see that the outgoing radiation 
is very similar to the incoming one. Fig. |^c. This case is similar to the ones in Fig. ^ and 
One can easily recover the "classical information", (the word "Hi"), from the outgoing 
radiation in Fig. Also the quantum radiation is quite symmetric. We can see from Fig. 
^ that the classical and quantum radiation are correlated. Actually one can recover the 
"classical information" from the quantum radiation alone. Each classical pulse results in a 
sharp decrease ("jump") of the quantum radiation. We see that for a pulse of a relative 
height 1 the jump is twice as large as the jump for a pulse of a relative height 0.5. Therefore 
by looking at the different jumps in the quantum radiation one can recover the "classical 
information" . 

In Fig. ^ the mass is about 80% of the critical mass, and the outgoing radiation is much 
more distorted than in Fig. As in the cases shown in Fig. and it is difficult to 
recover the "classical information" of Fig. ^c from the profile shown in Fig. 



B. Quantum correlation function 



As discussed earlier in section II, we choose the quantum state to be the vacuum on 53^, 
so that the quantum contribution to the energy-momentum tensor on is always zero, 
{Tyy){v) = 0. By construction (with the reflecting boundary conditions and no black hole 
formation), the evolution of the quantum state must be unitary. How can we see that the 
quantum radiation on is described by a pure state? The created quantum radiation 
reaching before the escaping classical matter is the beginning of the Hawking radiation 
which is almost thermal and by itself does not correspond to a pure state. A pure state 
can be recovered from the radiation on only if there are strong correlations between the 
early-time Hawking radiation and late-time negative-energy radiation reaching after the 
escaping outgoing classical matter. We demonstrate these correlations by calculating the 
correlation function 

c,.,M^,x') ^ {fl,{x)fl,Ax')) - (r/,(x))(r/,,,(x')). (45) 

The correlations between different points u and u' on are given by Cuu,u'u'iu,u'), where 
u is the asymptotically fiat null coordinate on 53^. For the reflecting boundary conditions, 
we have a closed form expression^ 

Cuu,u'u'{u,u) = — — — 4 , (46) 



Hn Ref. |T|], the factor Nh^ did not appear explicitly in Eq. (6) of Ref. and the values of 



[u, u 



quoted there are in units of Nh . 
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where vb{u) is the boundary curve in the {u,v) coordinates defined in Eq. (|3TD. Using the 
numerical solution for q{z), we numerically find the functions vb{u) = ln{q[z{u)]} and 
duVB{u), and from there can obtain Cuu,u'u'{u,u')- 

The denominator in Eq. (^) is a rapidly varying function of u and u'. By defining the 
relative correlation 

C( + ) _ Cuu,u'u'{u,u') 



C(-) Cu"u",u'u'{u",u') 



where vb{u) —vb{u') = vb{u') — vb{u") , (47) 



we get a measure for the correlations that varies more slowly as a function of u and u' and 
also is not proportional to h, which we take to be very small in the large N limit. This 
relative correlation is plotted as a function of u, with u' held fixed. The numerical results 



for the relative correlations (|47|) and the outgoing quantum radiation shown in Fig. ^a, when 
we fix Am' = 0, are shown in Fig. |^. Since for u < the quantum radiation in Fig. ^ is 
exponentially small, the correlation function C(— ), which describes the correlations between 
the radiation at u = and the radiation at m < 0, is approximately the vacuum correlation 
function (there are almost no created particles in this region). On the other hand, the 
correlation function C(+) describes the correlations between the radiation at m = and 
the radiation at m > 0. Since most of the radiation is in the region -u > 0, we can regard 
C(— ) as a reference function relative to which C(+) is being measured. If C(+)/C(— ) < 1, 
then the correlations described by C(+) are weaker than the vacuum correlations, while if 
C(+)/C(— ) > 1, then the correlations are stronger. 

We see from Fig. |^ that when u is in the region of positive-energy radiation (see also 
Fig. ^a), we have C(+)/C(— ) < 1. This is expected since the positive-energy radiation is 
the precursor of the uncorrelated thermal Hawking radiation. On the other hand, for u in 
the region of negative-energy radiation we have C{+)/C{—) > 1. Namely, the correlations 
between a point in the region of positive-energy radiation, i.e., at m = 0, and a point in 
the region of negative-energy radiation are stronger than the vacuum correlations. It is 
just when the energy of the quantum radiation becomes negative that C(+)/C(— ) becomes 
greater then one. Hence, the negative-energy radiation is not only necessary for energy 
conservation, but is also instrumental in recovering the correlations present in the final pure 
state. Before we can discuss what happens to information when M is very near the critical 
mass Mcr, we must study the properties of the near-critical spacetimes. 



V. NEARLY CRITICAL SPACETIMES 

A. Approaching the critical solutions 

The refiecting boundary conditions can only be imposed on those sections of the boundary 
curve that are timelike. Along these sections, the boundary curve can be written in the form 
), where well-defined function, and one can impose the condition 

(p!3|). Since the boundary curve is dynamical, the initial data determines its nature through 



the evolution equations. If the boundary is everywhere timelike, then the general solution we 
found earlier is regular, energy- conserving, and unitary. Such a solution is called a subcritical 
(i.e., a non-black hole) solution. If on the other hand the boundary becomes spacelike in 
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some regions of spacetime, then the spacetime develops a spacehke singularity [|T0[, which is 



initially hidden behind an apparent horizon. Such a solution is called a supercritical solution 
and describes an evaporating black hole [|I^. One can therefore distinguish between two 
regions in the space of solutions or the space of initial data: black hole and non-black hole 
solutions. Furthermore, one can find continuously varying parameters pi, specifying the 
initial data, such that all solutions with pi < p* are subcritical, while all solutions with 
Pi > p* are supercritical. The solutions on the boundary separating the two regions (i.e., 
the ones with pi = p*) are called the critical solutions. As a result of the non-linearity of 
the system, the physical properties of the solutions need not be continuous functions of the 
parameters at the critical values p*. 

We would like to study the transition between the subcritical and supercritical regions in 
our semiclassical theory. Using our numerical integration we can easily determine whether 
the solution is subcritical or supercritical: As long as the boundary curve is timelike we can 
write the boundary equation, = x~^{x~), in the inverse form ), where is the 

inverse function of x^. Since the boundary is smooth, the quantity dx]^/dx^ is continuous. 
For a timelike boundary dx~j^/dx~^ is positive, for a spacelike boundary dx'^jdx'^ is negative, 
and for a null boundary dx~^/dx~^ = 0. So at the point where the initially timelike boundary 
becomes null, just before becoming spacelike, the derivative dq/dz in Eq. ( P^D (which is 
the inverse of dx^/dx^) diverges and the numerical integration terminates. Therefore, all 
solutions for which the integration terminates at a finite value of are supercritical, while 
the others are subcritical. 

For the profiles in Fig. ^ the continuous dimensionless parameters in solution space 
are: M/A and Ae in case (a), M/X, Ae and Xvq in case (b), and M/A (the U are not 
continuous parameters) in case (c). We approach the critical solutions by varying one of 
the parameters, while keeping the others fixed. First we take m = M/X to be the free 
parameter. In case (a) we fix Ae = 5, and in case (b) we fix Ae = 0.1, Afo = 10. We find 
that in all the above cases there exists a critical value m* such that if m < m* the solutions 
are subcritical while if m > m* the solutions are supercritical. The value of m* depends 
on the values of the fixed parameters. To study the behavior of the solutions just below 
criticality, we numerically integrate the boundary equation for the cases < Am/m* << 1, 
where Am = m* — m. We take k = 10~^ in this section, in order to be able to easily probe 
the regime. Am < n. (The values of m near m* are of order 1 and thus are still large 
with respect to k.) Then we calculate the outgoing radiation on S^. In Fig. ^ we show the 
radiation on for the three cases (a), (b) and (c) of Fig. |^, where Am/m* ~ 10~^. We 
shift the M-coordinate to u — u^^, where u^^ corresponds to the last classical reflected null ray, 
(i.e., = X^^ ln[— Ax^(x2") — P+/A]). In Fig. ^ we show the total stress tensor, T/^, in units 
of A^, describing the outgoing radiation on S^J, and in Fig. ||b we show only the quantum 
part. The quantum radiation is scaled by 4/(fi:A^), such that {T^u) = 1 corresponds to the 
constant Hawking radiation from a 2D black hole ||19|| . 

We see from Fig. ^ that while the early-time radiation strongly depends on the details 
of the classical infalling matter, the late-time radiation in the region u> u^ — X~^ is almost 
the same for all the cases (a), (b) and (c). In the region u > — X^^ the classical part, 
T^^, increases rapidly to large values of the order of lO^A^, and then sharply decreases to 
zero. The quantum part, {Tuu), approaches the constant Hawking radiation, (T„„) = kA^/4, 
just before u = u^. It then decreases rapidly to negative values of the order of — lO^A^, and 
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finally increases rapidly to zero. This late-time behavior appears to be independent of the 
profile of the infalling matter, as long as < Am/m* << 1. Similarly, if we approach a 
critical solution by varying any one of the other parameters p, while holding the rest fixed, 
we find the same late-time behavior as long as < {p* — p)/p* -C 1. Let T^ax and T^m 
be the maximum and minimum values of the total radiated energy flux T/^ on For all 
the cases that we studied we find numerically that as {p* — p)/p* —>■ 0^, the value of T^m 
approaches —Tmax and Tmax tends to infinity. In the next subsection we analyze this critical 
behavior. 



B. Breakdown of the semiclassical approximation ? 



To understand the above late-time behavior analytically, we write T„^(m) and {Tuu{u)) 
explicitly, using (|1^ and (0), 




(48) 



and 



{Tuu{u)) 



1 ( dx \^ ( dx^{x ) \ ^ 



(Axg)2 \ du 



dx' 



(49) 



T^+, to be everywhere 



We take the classical stress-tensor describing the incoming matter 
regular, i.e., finite and smooth. Therefore the non-trivial contributions to (|^) and (^91) come 
from the "redshift factor", {dx~/du)^, or from the "blueshift factor", {dx^/dx~Y. Consider 
first the redshift factor. From the coordinate transformation (|3TD we get {dx~ /duf' = 
exp(— 2Au). If we ignore the back-reaction, i.e., take k = as explained in the boundary 
curve of the critical solution first becomes null as u —>■ oo, and the redshift factor approaches 
zero. However, when we include the back-reaction, the redshift factor at the point where 
the boundary becomes null is finite, as can be seen as follows. Let (x^,x~) be the point 
at which the boundary curve becomes null. At that point an apparent horizon is formed, 
since when the boundary curve becomes space-like (and describes a black hole singularity) 
it is surrounded by an apparent horizon |lO|. The equation for the apparent horizon is 
de~'^'^/dx~^ = [0, and using (^ we get 



Xx-^ix^) = P4x^)/X + 



4Ax+ 



(50) 



where the apparent horizon curve is x 
and using ( pT]) and (|50|), we find at (x 



^ah(^^)- Defining AP+ = P+(x"'" — > oo) — P+(x+) 



C 5 



X~ 



' dx 
du 



K 



4Ax: 



(51) 



Note that the apparent horizon is null at {x'^,x^ ). This can be seen by showing that the 
derivative of (^Dj) with respect to a;"*" vanishes at the point (x+, x~) on the boundary X = 0, 
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du - \ 4Ax+ / \ x+ / V 4Ax+ 



where X is given by Eq. ([21|). From the definition Pj^{x^) = J dx~^T'^_^_{x^) and the fact 
that dx~^/dx^ = at {x'^,x~), we get 

It follows that T^l (f (xjt)) is of order kA^, which is the same as the quantum contribution. 
Therefore, must be close to the cut-off point of the classical matter distribution (the 
location of is shown in Fig. assuming that the profile T^l (v) has no anomalously long 
tail. Expanding AP+ ~ {dP^/dx~^)Ax~^ = T^'_^(x"'")Ax"'", where Ax"*" = X2' — x^ -C x^, and 
using (|50D-(|52D we finally get 

( dx^\^ ( K \^ ( ^ 2Ax+\ / K \^ 

(53) 

Thus the redshift factor (|53D is finite at (x+,x~), and is zero only if we neglect the back- 
reaction (k = 0). 

Next consider the blueshift factor, {dx^/ dx~Y . Obviously when the boundary curve 
becomes null the blueshift factor diverges. This divergence, together with the finiteness of 
the redshift factor and is the reason for the divergences in (^Sf ) and p9|) at Uc = u{x~). 
But what about the total stress-tensor, T/^ = T^l + (T„u)? From (^Sp, (^), and (0) we see 
that at (x^, x~) the diverging terms in T^^ and (T„u) cancel each other, leaving a finite T/^: 

TlM = ^ . (54) 

So why do we get divergences in the total stress tensor, shown in Fig. |a? To understand 
this behavior, consider a point (xq , Xq = x~^{xq)) near (xjt, x~) on the boundary of a critical 
solution. Defining 6 = x'^ — Xq and expanding (^8]), (^9]) about (x^,x~) to leading orders 
in S, dropping the finite term (pi), we find 



4Axj 




' (55) 



{a6 + M2)2 



where a = [d'^x]^/d{x^Y]c, b = [d^x]^/d{x^Y]c, and uq = u{xq). The point (x^,x~) cannot 
be a local maximum of the boundary curve, since the boundary curve of the critical solution 
never becomes space-like. It also cannot be a local minimum since by construction it is the 
first (and only) point at which the boundary is null. Therefore the point (x^,x~) must be 
an infiection point, i.e., a = 0. Also, for a general critical solution we have {^dT'^^ / dx^^ ^ < 

— k/[2(x^)3], which is necessary for the apparent horizon in (|50D to become timelike beyond 
the point (x^,x~). We therefore get from (|55D that T/^(mo) — A5~^ with A > 0. For 
points on the boundary to the past of x^, T/^ ~ and for points to the future of x^, 

~ — A|(5|~3 as \5\ — > 0. Thus for critical solutions the radiation fiux diverges in opposite 
ways on both sides as x^ ■ For almost critical solutions one gets the nearly divergent 
results shown in Fig. |^. 

The fact that T/^ diverges does not necessarily signal a breakdown of the semiclassical 
approximation in which the metric and dilaton fields are treated as classical dynamical fields. 
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From conservation of energy we know that the total amount of outgoing radiation (the sum 
of positive and negative parts), is finite and equals the total amount of incoming radiation. 
Moreover, the total amount of positive-energy radiation alone, E^, in the region between 
x~ and Xq on 9^^, is proportional to Tmax^x" , where Ax~ = x~ — Xq . From x+ — Xq = S, 
we find that Ax~ = bS"^. It then follows from Eq. ( |55D that is finite. The same holds 
for the total amount of negative-energy radiation. Thus the divergence in the density, T/^, 
can be viewed as describing a "shock wave" or "thunderpop" [^. Although a shock wave 
does not necessarily imply a breakdown of the semiclassical approximation, an examination 
of the fluctuations in T^^ provides strong evidence for such a breakdown. Let us calculate 



the fluctuations in T^^^. As an estimate for the fluctuations one can use the quantity |]33|j34 



^ ^ujj.'u'i,-^^ 2; ) 



(T^/.(x)f/,,,(x')) - (T/,(x))(f/,,,(a:')) 



c, 



(x, x') 



{fUx)){ffix')) 



(56) 



Far below criticality, with x and x' sufficiently apart, the numerator in Eq. (EBf) will be less 



than or of the order of A^^ (see Eq. (^61) ). On the other hand, the denominator is of the 
order of = N'^h'^. Hence, in this case the fluctuations (|56|) will be of the order of 1/A^ even 
in the region of negative energy. In the large A^ limit these fluctuations are small and the 
semiclassical approximation is valid. However, for critical or nearly critical solutions this is 
no longer the case. Let us estimate A'uuu'u'iuyu') for a critical or nearly critical solution, 
with u and u' sufficiently far apart that vb{u) — vb{u') ~ which sets a natural length 
scale. Let u' approaches Uc, and take u to be in the region of negative energy such that 
vb{u) — vb{u') ~ A~^, which is consistent with the quantum inequality ([371) . Then using 
and (pBf), we get for u' —>■ 



A', 



1 / SvbW) 



(57) 



For any given A^, no matter how large, while keeping u fixed, we can find u' sufficiently 
close to Uc such that the blueshift factor, {OvB/du)"^, in (|57D is dominant over A^, and the 
fluctuations (|57|) are large. We regard this anomalous behavior as implying a breakdown of 
the semiclassical approximation at criticality. Unlike in the 4D case, where the temperature 
of the black hole formed just above criticality approaches infinity, in the 2D case the black 
hole temperature is always a finite constant, Tbh = A/27r, but nevertheless the semiclassical 
approximation (even in the large A^ limit) breaks down at criticality. Also, for the critical 
solution the quantity p' = dx']^/dx~ diverges at x~, and dropping the last term on the r.h.s. 



of Eq. ([T5|) , i.e., the moving mirror term, becomes problematic even in the large A^ limit. 



VI. HIGH DENSITIES OF INFORMATION 



A. Quantum limitations on the negative-energy radiation 



Let us calculate the relative-correlation function (H^ for nearly critical solutions. As in 
Fig. we take the infalling matter to be the one shown in Fig. with Ae = 5, but while 
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in Fig. 1^ we take Am/m ~ 0.9, in this section we take Am/m ~ 10~^. The results for 
the relative correlations are shown in Fig. 03- We see that in the region of the Hawking 
radiation (0 < Am < 10.5) the relative correlations are very small. Then in the region of the 
brief burst of negative energy, the relative correlations increase sharply to extremely high 
values of the order of 10^^, and sharply decrease to one. It is this sharp increase in the 
relative correlations that insure a unitary evolution. If it were not for the negative-energy 
burst the correlations would be lost and the final state would be a mixed state. 

The information encoded in these correlations is very dense as the burst of negative 
energy is very localized (see Fig. ^a). To find an upper bound on the duration of the 
negative-energy burst, Au, we use the quantum inequality (^). To estimate the value of 
|T/„|, we use the analytic expressions (|48|) , (|49|) for a nearly critical solution. The quantity 
|£'neg| is approximately of the order of 

l^negl ~ \{TlU^\Au ^ (T/J„,,.An . (58) 

From (^Sp and ( PD| ) one can see that {Tl^)max corresponds to the maximum of dx'^/dx'. 
Let ( 

^m?-^m) point on the boundary curve for which T^^ is maximum. Using the 

boundary equation X{x^,x]j) = Xb and the solution (^) we find that at the maximum of 
dx% I d: 



X 



- ^ ■ (59) 

From (^HD, (^91), (|53[), and ([59| ) we find that for a nearly critical solution (for which 
(5x^/(9x"'")mm is very small) 

\ 171/ \ / mm 

From (|58|), (pOD and the quantum inequality (0), we get 

Au<^-^(^] . (61) 
K, \ dx+ I . 

\ / mm 

Since x^ < X2 is finite, as one approaches the critical solution {dx^/ dx^)min approaches zero 
and so does Am. If the amount of information. A/, carried by this burst of negative energy 
is finite, then as we approach the critical solution the information density, / = A J/ Aw, 
diverges. In the next subsection we show that this is indeed the case. 



B. Entropy and Information 

The results regarding the quantum part of the energy-momentum tensor discussed in 
the previous sections may plausibly be interpreted as arising from the creation of particle- 



antiparticle pairs ||35[. The particles reach infinity, Sj, and give rise to the positive-energy 
(Hawking) radiation, while the antiparticles carrying negative energy are reflected from the 
boundary and give rise to the negative-energy radiation on SJj. If it were not for the negative- 
energy burst of radiation, the correlations between the particles and antiparticles, shown 
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explicitly in Fig. 0, would be lost and the final state would be a mixed state with non-zero 
entropy. The entropy of this mixed state can be found using the fact that the spectrum of the 
outgoing positive-energy quantum radiation is almost thermal, as can be seen by calculating 
the corresponding Bogolubov coefficients [^. The entropy of this thermal radiation is the 



dimensionless Boltzmann entropy [21 



SboI = 2h~'^ ^/S{u}du , (62) 



where S is the energy density of the thermal quantum radiation, S = {Tuu), and mq is the 
value of u at which the thermal radiation ends. From Fig. ^b we see that as we approach 
the critical solution, the thermal quantum radiation is almost independent of the specific 
profile of the infalling matter. Therefore to get an estimate of (^) one can calculate it for 
the case of very localized infalling matter. In that case one finds analytically that S is given 
by 

^ = ^{'- (l + XAe^^r ) ' ^^^^ 

where A = M/{X^Xq), and = Xq is the null trajectory of the localized (shock-wave) 
infalling matter. Then the Boltzmann entropy of the thermal radiation is 

S^,^%m], (64) 



6 V Ak 

where we assume that ln(4M/A/t) ^ 1. The entropy (|6^) is by definition the amount of 
information that is lost by ignoring correlations between the thermal radiation and the later 
burst of negative energy, Jiost = 5'boi- This information cannot be recovered from the thermal 
radiation before the arrival of the negative-energy burst. After the burst of negative energy 
all the correlations between the particles and antiparticles are restored, and we get a pure 
state of zero entropy. This is expected since we impose reflecting boundary conditions and 
our spacetime has a trivial topology. Therefore, the information A/ that is gained during 
the arrival of the burst of negative energy equals /lost, 

AI = /lost = 5bo1 ^ ^ In . (65) 

As we approach the critical solution, AI remains finite. Using (|6T|) and (65) we get 



• AI iVA, /4M\ K fdxt\ 

A^^ 4^^Hl^J(3<Fl^j • ^ ^ 

v my \ / max 

As the critical solution is approached, / — > oo. 

Even though our space is one dimensional, the result (|66D suggests that the way informa- 
tion is transferred by the negative-energy burst is different from that by which information 
is transferred in linear channels. The theoretical upper-bound on the bulk of information 
flow in linear channels is 
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/linear <^ln(iV), (67) 

where the "message" is transferred with positive energy E. Suppose that this bound could 



be extended to negative energies, £" < 0, by replacing E in (67) by \E\. From Eq. (|6^) , we 
find 

\ / A T\/r \ 

(68) 

Then the theoretical upper bound ( |67D for our system would be 

NX , /4M\ , ,,,, 

- In(A^) . (69) 





fAM\ 


l^l^^ln 






\ k\ ) 



/linear — 1^ 



967r \ k\ J 



For any given sufficientely large value of N one can find a range of solutions just below 
criticality for which the semiclassical approximation remains valid (i.e., the fiuctuations 
(^) are small) while the bound (|69D on /unear is exceeded by I in (|66|) . This is achieved by 
requiring that a/ZV >> {dvB/du)^^,^ > ln(iV)/87r. 

It may be better to view the information associated with the negative energy in terms 
of informationstorage rather than information transfer. The brief burst of negative energy 
can be considered as a very localized configuration containing a finite amount of information 
that is traveling in space. As one approaches criticality the density of stored information 
becomes unbounded. This seems to be in agreement with other considerations regarding the 



unboundedness of information storage densities in quantum systems 38 



VII. CONCLUSIONS 

In this work we present and study in detail a theory of semiclassical 2D dilaton gravity 
with reflecting (conformal) boundary conditions. This theory shares many of the features of 
spherically symmetric semiclassical gravity. In particular, massless scalar fields can collapse 
to form a black hole that evaporates. If the energy and energy density of the initial configu- 
ration of the scalar fields are below certain critical values, i.e., the subcritical or weak initial 
data case, then the scalar fields do not form a black hole, but instead are reflected from the 
boundary and escape to infinity. On the other hand, if the energy and energy density are 
above the critical values, i.e., the supercritical or strong initial data case, then the scalar 
fields collapse to form a black hole. This created black hole evaporates by emitting Hawking 



radiation. In a previous work |jTO[ we found that such evaporation leads to an end-state 
geometry similar to that of a 4D semi-infinite throat. Here we study subcritical solutions, 
especially those near criticality, for general smooth initial data. 

By combining analytical and numerical techniques we investigate the detailed structure of 
the classical and quantum one-loop contributions to the outgoing radiation reaching asymp- 
totic future null infinity, For the subcritical solutions, we find that before the refiected 
massless classical matter fields reach 9^, positive-energy quantum radiation is observed at 

and continues as the classical matter reaches 53^. This positive-energy quantum radiation 
is followed by a fiux of negative-energy quantum radiation. 
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For profiles of incoming classical matter that encode information in bits consisting of 
pulses of two different amplitudes, we calculate how the encoded information appears in the 
outgoing classical profile on S^J. For the quantum radiation, the situation is more involved. 
The outgoing positive-energy quantum radiation by itself does not describe a pure state, but 
is strongly correlated with the negative-energy quantum radiation reaching at a later 
time. This negative-energy quantum radiation not only insures conservation of energy but 
also restores the correlations necessary for a pure final quantum state. 

The study of solutions just below and at criticality gives insight into the information 
puzzle in a unitary framework. As one approaches the critical solution the classical outgoing 
radiation becomes very distorted at late times. Part of that classical radiation reaches 9j 
with a time-delay and has the form of an extremely dense brief pulse. As the critical solution 
is approached the late-time classical energy density becomes highly distorted and ultimately 
diverges, making it impossible to recover the complete classical information. Regarding 
the quantum radiation, as one approaches the critical solution the early-time radiation 
is indistinguishable from thermal Hawking radiation from a black hole, while the energy 
density of the late-time negative-energy radiation diverges, making recovery of the quantum 
correlations impossible. Although the late-time energy density of the outgoing radiation is 
infinite in the critical case, the total amount of energy is finite and conserved. Nevertheless, 
the semiclassical approximation breaks down because the fluctuations in T^^ become very 
large. 

The above analysis shows that the black hole phase transition (subcritical — > supercrit- 
ical) and the information puzzle are intimately related. At the critical solution there is an 
apparent breakdown of predictability. However, this breakdown of predictability is related 
to the breakdown of the semiclassical approximation. The same divergence in energy density 
that makes recovery of information and quantum correlations impossible evidently makes 
the semiclassical approximation invalid due to large fluctuations in T^^^. Wc would like to 
stress that even though we take the large limit to enforce the validity of the semiclassical 
approximation for subcritical solutions, the semiclassical approximation nevertheless breaks 
down at criticality. This may suggest that the breakdown of the semiclassical approxima- 
tion at the onset of black hole formation is a fundamental result independent of the speciflc 
model to be studied. This is supported by the fact that the crucial features discussed in 
Sec. V B, i.e., the flnite redshift factor and inflnite blueshift factor, seem to be independent 
of the explicit semiclassical model. Thus, based on the semiclassical approximation alone, 
it may not be possible to trace the fate of information and correlations at the onset of black 
hole formation. The full quantum theory appears crucial to understanding if the overall 
evolution can be unitary at and after the onset of black hole formation. Nevertheless, our 
results do show that as criticality is approached via subcritical solutions, information and 
correlations become essentially irrecoverable while unitary evolution is preserved. 
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FIGURES 



FIG. 1. Penrose diagram of the Linear Dilaton solution. The heavy curve is the boundary 
curve separating the regions of weak and strong couphng. 

FIG. 2. Penrose diagram of a typical subcritical solution. The infalling matter with support 
on the interval < < .t^ on 3^, is reflected from the boundary to become outgoing toward 
future null infinity, Here z = \x~ . 

FIG. 3. The incoming stress tensor, Ty^{v), for three different families of initial data, with (a) 
corresponding to Eq. (32), (b) corresponding to Eq. (38) and (c) corresponding to Eq. (40). 

FIG. 4. The stress tensor T^u in units of describing the outgoing radiation on G+. The solid 
curves describe the total outgoing radiation, T^^, while the gray curves describe only the quantum 
part, {Tuu), scaled by 2/k. In (a) the mass is about 6% and in (b) it is about 65% of the critical 
mass for the infalling matter shown in Fig. 3(a). In (c) the mass is about 3% and in (d) it is about 
65% of the critical mass for the case shown in Fig. 3(b). 

FIG. 5. The stress tensor T„„ describing the outgoing radiation for different values of Ae in the 
case of the infalling matter shown in Fig. 3(a) with M/A = 0.1. See the text for more details. 

FIG. 6. The stress tensor describing the outgoing radiation for the infalling matter shown 
in Fig. 3(c). The solid curves describe the total outgoing radiation, while the gray curves describe 
only the quantum part. In (a) the mass is about 4% of the critical mass, while in (b) it is about 
80% of the critical mass. 

FIG. 7. The relative correlation function of the outgoing quantum radiation on In (a) the 
mass is about 6% of the critical mass, while in (b) it is 99.999% of the critical mass. 

FIG. 8. The stress tensor T^^ describing the outgoing radiation for nearly critical solutions 
with M being 99.999% of the critical mass. In (a) we plot the total outgoing radiation and in (b) 
only the quantum part is plotted. The dashed, dotted and gray curves correspond to the infalling 
matter shown in Figs. 3(a), 3(b) and 3(c), respectively. 
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FIGURE 1 




Fig. 1: Penrose diagram of the Linear Dilaton solution. The heavy curve is the boundary curve separating the regions of 
weak and stron couphng. 



FIGURE 2 




FIGURE 3 




Fig. 3: The incoming stress tensor, Tyy, for three different famihes of intial data, with (a) corresponding to Eq. (32), (b) 
corresponding to Eq. (38), and (c) corresponding to Eq. (40). 



FIGURE 4 




Fig. 4: The stress tensor, Ty^y^ in units of describing the outgoing radiation on The sohd curves describe the 

total outgoing radiation, T^^, while the gray curves describe only the quantum part, {Tyy^ ^ scaled by 2//€. In (a) the mass 
is about 6% and in (b) it is about 65% of the critical mass for the infalling matter shown in Fig. 3(a). In (c) the mass is about 
3% and in (d) it is about 65% of the critical mass for the case shown in Fig. 3(b). 



FIGURE 5 




Fig. 5: The stress tensor Ty^y^ describing the outgoing radiation for different values of A 6 in the case of the infalhng 
matter shown in Fig. 3(a) with M/A = 0.1. See the text for more details. 



FIGURE 6 




Fig. 6: The stress tensor describing the outgoing radiation for the infaUing matter shown in Fig. 3(c). The sohd 

curves describe the total outgoing radiation, while the gray curves describe only the quantum part. In (a) the mass is about 
4% of the critical mass, while in (b) it is about 80% of the critical mass. 
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Fig. 7: The relative correlation function of the outgoing radiation on In (a) the mass is about 6% of the critical 

mass and in (b) it is 99.999% of the critical mass. 
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Fig. 8: The stress tensor with M being 99.999% of the critical mass. In (a) we plot the total outgoing radiation 

and in (b) only the quantum part is plotted. The dashed, dotted and gray curves correspond to the infalling matter shown in 
Figs. 3(a), 3(b), and 3(c) respectively. 



